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Abstract: We recast the study of a closed string gas in a toroidal container in 
the physical situation in which the single string density of states is independent 
of the volume because energy density is very high. This includes the gas for the 
well known Brandenberger-Vafa cosmological scenario. We describe the gas in the 
grand canonical and microcanonical ensembles. In the microcanonical description, 
we find a result that clearly confronts the Brandenberger-Vafa calculation to get the 
specific heat of the system. The important point is that we use the same approach 
to the problem but a different regularization. By the way, we show that, in the 
complex temperature formalism, at the Hagedorn singularity, the analytic structure 
obtained from the so-called F-representation of the free energy coincides with the 
one computed using the S-representation. 
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1. Introduction 

Treating a gas of free closed superstrings when all the spatial dimensions are closed 
and the system is kept in a thermal bath leads us to the well known conclusion that 
the Helmholtz free energy diverges as one approaches the Hagedorn temperature 
from the low temperature regime. This way, as the energy U also diverges, one can 
conclude that, in a fixed temperature description, the Hagedorn temperature is a 
maximum one for the system. Let us recast in the next section what the details 
for the description in the macrocanonical ensemble with /i — are. In particular, 
we will remind to the reader that energy fluctuations give us a measure about the 
very possibility of a fixed temperature description of the system. In section 3 we 
will compare our results with the thermodynamics that stems from the generalized 
ensemble description of an analogous extensive system. After all this, we will have 
then got some results to get into the treatment of the system at fixed energy in 
section 4. There, in a first subsection, we will remember and critically recast the well 
known computation by Brandenberger and Vafa [1]. Next, we will present another 
computation that will reinforce the conclusion that the specific heat, as a function 
of energy, is divergent. A final subsection will be devoted to explain whether volume 
dependent corrections can change the picture. Finally, section 5 will present a few 
comments and serve as a reminder of the main results. 



2. The macrocanonical description of closed strings at finite 
size 

In a macrocanonical description at null chemical potential, the grand canonical parti- 
tion function can be equated with the canonical partition function at a given number 
of strings N*. This number maximizes the canonical partition function Z([3, R, N) 
(i.e., [d Z/d N] N=N , = 0). When the fluctuations in the number of strings are 
small, this maximum coincides with the averaged number of strings, N, as com- 
puted in the grand canonical ensemble. To be more concrete: (/3, R, fi = 0) = 
YHv=o ^ (P, R, N) « Z (/3, R, N* R)) and Z finally results a function of (3 and R 
only. This Z is what we call the partition function for the system of free strings. 
Minus the logarithm of the partition function divided by (3 is what we call the free 
energy of the system and it is only a function of R and /3 and not of the number of 
strings (see [2]). This will exactly correspond to the thermodynamical free energy 
whenever a thermodynamical limit can be defined. The computation for the gas of 
free superstrings gives 
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This is the free energy in the so called ^-representation. It is the result one gets 
when computing the Helmholtz free energy in the light-cone gauge or by summing up 
over the field content of the string (analog model) (see [3]). In the conclusions we will 
comment more on this point and its relation to the F-representation, that coincides 
with the computation of the free energy as a vacuum energy for the Euclidean theory 
with Euclidean time of length /3 including winding modes around it. 

Now, let us introduce an ultraviolet dimensionless cutoff e in r 2 . This automat- 
ically produces the splitting of the free energy as F (j3) = F l {(3) + F h (j3) where F l 
is got 1 by integrating r 2 from e to +oo and F h by integrating the same integrand 
over r 2 from to e. 

The integral over T\ simply represents the left-right level matching condition, 
i.e., a Kronecker delta of the form #jv-jv+m-n,o- Here, TV and N stands for the right 
and left oscillator numbers. 

Supposing e C 1 has various implications. The Jacobi # 2 function, having the 
/3 dependence, can then be approximated by the leading terms (two in fact) in the 
series expansion representing it. This is physically equivalent to taking the classical 
( Maxwell- Boltzmann) statistics approximation. One can also use the fact that r 2 < e 
to approximate with arbitrary precision the contribution from integrating over Ti, 
i.e., implementing the left-right level matching condition. 

To compute first the integral over t±, it is very useful to note that, when r 2 
is small, the main contribution to the integrand will come from a neighborhood of 



Furthermore, one can choose e small enough so as to have that ea' / R 2 <C 1 and, 
simultaneously, eR 2 /a' <C 1. It is enough to choose e to fulfill the most stringent 
criterion. In fact, one criterion converts into the other by T-duality (R — > a'/R). If 
R > \fa! ', one can choose e a'/R 2 in order to accomplish both criteria. When 
R < \fa' we have that e <^ R 2 /a' suffices to satisfy both relations. We will see 
how these criteria, that define the physical system we are treating and the kind of 

1 F l , as it is written in this proper time representation, shows a divergence when T2 — > +00 for 
the massless excitations of the string when momentum and winding numbers vanish. This is an 
artifact of this representation that results from the second quantization of the vacuum state then 
producing a divergence as lnr2 when T2 goes to infinity. This divergence should not be present with 
finite volume or, more precisely, when the momenta are not dense, and can be subtracted without 
affecting the results of our work. 
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thermodynamical limit we are taking, can be expressed in terms of energy density. 
This allows us to approximate the sum over m and n (i.e., over winding and mo- 
mentum numbers) by a multiple integral over the whole M 18 on the real variables 
mi, ■■■, m 9,ni, ...,n 9 . This is a valid approximation as given by the Euler-Maclaurin 
formula (see [4]). The multiple integral over windings and momenta shows us that, 
for the contribution coming from the ultraviolet degrees of freedom to F ((3) (i.e., 
the contribution encoded in F h ), no dependence on R, and then on the volume 
V = (2nR) 9 , will survive. This is so because the change of variables of unit Jacobian 
m — ► m\fa'/R, n — ► nR/yfa' makes R to disappear in the multiple integration. 
Everything happens as putting R = y/a' that simply shows the fact that a situation 
in which R is in a neighborhood of \fa' also suffices to compute the sum by an 
integral. This will be linked to the cosmological situation in which all the spatial 
dimensions are of the order of the selfdual length and the system is, in some sense, 
small. 

With all this together one can finally write 



+1/2 
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Where the modular properties of the transverse partition function have been used. 
We remark again, that the main contribution in the r 2 — > limit for the T\ integral 
comes from the neighborhood of T\ — 0. The coefficients Oj are computable numbers 
(see the Appendix). In particular, clq = 1 

F h ((3) can then be approximated, for 2-n ((3 2 - (3^) <C e(3^ , as 2 

_^.r(,^)|.(^)- ( ,3) 

Ph = 7rV8 a' is the inverse Hagedorn temperature for closed superstrings type IIA 
and IIB (both have the same free energy because are indistinguishable at finite T) 
and (3 — Ph > 0. The b n are coefficients that can be directly connected to the a,; for 
example, bo = 1/(3h, and, in general, the b n oc 1//3h are independent of e computable 
numbers. The important point now is that, as discussed in the Appendix, 

n=0 \ PH J n=0 PH P 



2 The well known relation T [a + 1, x] = aT [a, x] + x a e x is very useful to this purpose. 



- 3 - 



As it is implied by the behavior of T[0, z] for big z, F h goes exponentially to zero 
as (3 — (3h grows. Because there is no dependence on the volume, it is not true that 
we can write F (f3) = —PV. This gives more importance and justifies the detailed 
way we have introduced F (f3) at the beginning of this section. It is now clear from 
the behavior of r[0, z\ when z goes to zero that the contribution of F h to the free 
energy when (3 ~ f3 h is much bigger than that of F l because \F h \ grows unbounded 
as (3 — > (3ft as long as F l gets the finite value F l (p H ). 

The concrete behavior around (3h can be made more explicit by using that f3 2 — 
{31 = 2(3 H ((3 - (3 H ) + (P~p H ) 2 ^ 2(3 h {(3 - (3 H ). One finally gets 
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This, by using inverse Laplace transformation, easily provides us the main ingredient 
to get the fixed energy description, Qi(E), i.e., the single string density of states. 

Sl*(E) = 9(E-A) ^ (2.6) 

The step function shows the utility of the dimensionless e parameter by imposing 
the condition E > A = 47r/(e/3#). This, when R > y/a', finally enforces E 3> 
A-kR 2 /(f3 H a') as the condition for the validity of (2.6) 3 . For the T-dual situation, 
one gets the dual condition. 

It is now an immediate task to get U {(3) around (3h 

rrh (as = 9 [P Fk (P)] = 1 -47r(/3-fa)/(6/3 g ) 1 (n 7) 

yf) ~ dp " P-Ph P-Ph k ' j 

Where we have taken that < P - p H < eP H /(Aix) = 1/A. 

Fluctuations in the macrocanonical energy are very useful for our study. They 
can be computed to give, near Ph 

[r'ct(T,v)Y 12 = 1 + o[{0 _ /3h)1] (28) 

So energy relative fluctuations are finite and not negligible. This means that we 

may expect the fixed energy ("micro") description and this macrocanonical picture 

not to be equivalent. 

3 It is important to remark that E is here the energy of the states which are accessible by one 
string. This condition can be compared to the one in [4] for the validity of converting the sum 
over winding and momenta into an integral in the direct calculation of Ti from its very definition 
(fix = dTi/dE); the condition is E > R/a'. It is clear that E > AnR 2 / (/3 H a') implies E > R/a' 
if R > Vrt. 
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The entropy can easily be computed as a function of (3 to give 



The fundamental relation in the entropic representation can now easily be ob- 
tained to be 

S h & (3 H U h + \nU h (2.10) 

In the entropic representation it is manifest that the non extensive term lnlf h 
is the one responsible for the positivity of the specific heat. When energy is high, 
the logarithm of the averaged energy is very small as compared to energy itself. If 
one sees the big fluctuations as giving the error of the energy variable to produce 
S(E = U), one perhaps should neglect the logarithmic term because is smaller than 
the error. 

The macrocanonical calculation gives us the number of strings as — f3F {(3) if 
Maxwell-Boltzmann statistics applies. Namely 

—/3F h {0) = N h ^ -In [((3 - (3 H ) A] ^lnU h (2.11) 

One can write the entropy in terms of the number of strings as a function of 
energy to get 

S h ^NiU^+e^ 11 ^ (2.12) 

This expression for the entropy can be compared to the one for a regular system 
for which the entropy scales with the number of objects as in the black body, for 
instance. In our gas we have that the entropy near (3h (that also gives high U) grows 
exponentially with the number of objects. This behavior in terms of the number 
of strings can also be compared to the open superstring case in which T H is a true 
maximum temperature [6]. For it, the entropy grows as 2N{U h ) + KN (U h )/V, 
with K a constant, and is a degree one homogeneous function of energy and volume 
(extensivity is a property of the gas of open superstrings in the infinite volume limit). 

It is then clear that the specific heat, as a function of the temperature, is positive 
for the high temperature phase, i.e., when T is near T#. The problem is that the 
order one energy fluctuations tell us that U h is a bad canonical average for the high 
energy of the non-isolated system. 



3. The generalized ensemble and extensivity 

The description at fixed temperature we have made is one very special. The reason 
is that our system does not depend on volume because this variable does not appear 
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in the description of the system which is also at fixed temperature and null chemical 
potential. 

A description through a generalized ensemble is one in which the system is char- 
acterized by intensive parameters. In a simple system they are pressure, temperature 
and chemical potential instead of volume, energy and the number of objects which 
are the corresponding extensive parameters. Since our string gas is one for which no 
volume dependence appears 4 , the ensemble we have named grand canonical is really 
a generalized ensemble. 

This ensemble does not appear thoroughly treated in regular textbooks on Sta- 
tistical Mechanics but can be found in [7]. The main point to take into account when 
reading this textbook is that the treatment of the description using this ensemble 
depends on the fact that extensivity is assumed for the system. The first notable 
fact about this ensemble, if extensivity is assumed, is that fluctuations in volume, 
energy and the number of objects must be big. The reason is that, in this picture, 
the system is characterized by pressure, temperature and ji and then volume, total 
energy and the number of objects can get any value with equal probability. It is 
worth to remark that this must be so when extensivity holds. 

Assuming extensivity in our problem would imply j3F (j3) = because Euler's 
relation (U — TS + PV — jiN = 0) would hold and we have /i = and no volume 
dependence. On the contrary, non extensivity allows a non vanishing free energy 
as computed previously, the relative energy fluctuations to be big (order one) and, 
simultaneously, the fluctuations in the number of strings to be small when energy 
is big enough. Indeed, the fluctuations in N(T) are small because, when Maxwell- 
Boltzmann statistics and the dilute gas approximation hold, they are given for any 
system by 

Van* _ \l{zd z ) 2 {zg) i _ i 

where q = —f3F (f3) is the single object partition function, a function of T and V in 
general and only of T in our problem. For the classical counting, q gives the number 
of objects. The general application of this result is not in contradiction with the 
fact that, in the generalized ensemble when extensivity holds, the fluctuations in the 
number of objects are big. When the system is extensive all the subsystems in the 
generalized ensemble with different object numbers are equally probable and then, 
to get the partition function, the sum over the number of objects does not run up to 
infinity (see [7] again). 

4 This is different from being a problem at zero pressure. In a system at zero pressure, the volume 
would be a function of temperature. In our case, there is no volume dependence at all and then the 
pressure vanishes. 
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On the other hand, the fact that for any extensive system at [i — and with no 
volume dependence one has U — TS = has immediate consequences. By putting 
T = dU/ dS, U — S dU / dS = can be understood as a differential equation that can 
be solved to give S = (3 U where (3 is a constant of integration that gives the constant 
inverse temperature of the system. In our problem this fio is Ph- From the point 
of view of Legendre transformations, this is the most simple and extreme case for 
which the transformation cannot be defined because the system has an infinite specific 
heat. Furthermore, if we now associate a density of states to the obtained entropy 
we have Q (E) = Kj3oe l3 ° E , where K is a positive dimensionless constant. Computing 
now back the partition function Z (f3) thorough Laplace transformation we obtain 5 
Z ((3) = K{3 J+°° dE e -^~M E = K(3 / ((3 - (3 Q ). So we find the surprising result 
that the free energy is not zero, but actually (3F (f3) = In [((3 — f3o) / (K(3 )]. The 
free energy must be negative and then has physical meaning when (3 is near and 
bigger than (3 so, after all, there is a maximum temperature T and there must also 
be an energy cutoff for the validity of Vl (E). Finally, one exactly gets the singular 
dominant term in the incomplete Gamma function in (2.5) identifying To = T#; i.e. 
the logarithmic contribution. This is a very simple example of non equivalence of 
ensembles so brightly explained in [8] and that relies upon the interplay between 
Legendre and Laplace transformations. It is very important to notice that the single 
object partition function — f3F (f3) we get actually coincides with the one in (2.11). 

Conversely, one can easily show that if the canonical energy is a function of T 
such that, at T = T , ((3 — f3o) U (f3) gives zero (or a constant) when (3 approaches 
Po then, since S(f3) = (3U((3) for a volume independent extensive system at \i = 0, 
one gets that S(f3) = (3 U((3) around (3 . But, in fact, we have already arrived at the 
fundamental thermodynamic relation S = (3q U that implies that the temperature is 
constant and equal to T . This contradicts the hypothesis assuming that the internal 
energy is a function of a variable (3. The final output is then that S = /3 U only makes 
sense in microcanonical thermodynamics and, more than this, extensivity holds for 
the microcanonical thermodynamics we get. 

It is certainly a notorious fact that, in this particular case, one can deduce the 
exponential growth of the density of states with energy from the hypotheses that 
extensivity holds, the system does not depend on volume and equilibrium is got at 
zero chemical potential. In fact, all the above reasoning simply tells us that £l(E), 
the density of states for the gas of strings, is a constant times e 13 " E when energy goes 
beyond a certain value, let's call it 0. However, this contradicts the computation in 
[1] in which a different non independent of the volume non extensive high energy 
entropy is found. This makes us suspect that there is something wrong in that 
computation. 

On the other hand, the grand canonical description of the closed string gas near 
5 We could have included a cutoff <\> to get Z (0) = K /3 e-<^-^7 (0 - O ). 
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Th and under the condition that windings and momenta are equivalent is really a 
generalized ensemble description but without assuming extensivity. Non extensivity 
is usually related to a kind of smallness of the system in size or number of objects, 
the presence of long range interactions or a critical behavior. We will dwell a bit 
more on this point in the final section. 

4. The fixed energy description of closed strings at finite size 

In order to study the fixed energy description of the gas of closed strings, an expres- 
sion for the multiple string density of states, Q (E), is needed. This has been done in 
the past in several ways. One could use, for example, a saddle point approximation 
to calculate the asymptotic expression for the density of states of the string gas, 
using the fact that f2 (E) = {exp (— f3F {(3))}. However, once the single string 
density of states is provided, Q (E) can be obtained using the convolution theorem 
[9]. This is the method used long time ago by R. Brandenberger and C. Vafa [1] (see 
also [10]), who calculated 6 



fli(E) can be expressed as sum of two terms Q[{E) and (E) where the superscripts 
/ and h refer to the low and high energy expressions. This introduces a cutoff A 
separating both regimes that can exactly be written as our cutoff An/^ePn) on the 
energy of one string. In fact we can obtain Qn(E) with a high degree of accuracy only 
from the convolutions of &i(E) as in the gas of open strings [6]. The convolution 
between Q[(E) and &i(E) to give a contribution to Vt^E) is negligible, as can 
clearly be seen in Fig. 1 where the comparison between the exact calculation 7 of 
u>2(E, t) = Qi(E — t) Vt\{t) and its approximation resulting from considering only the 
high energy part approximated by in (2.6) is shown at E = 56, R = Vc7 = 1. 002 
is a measure of equipartition of energy between the two strings. In fact, it is Fig. 1 
that lets us fix the cutoff A = 4r- 



4.1 The calculation of Brandenberger and Vafa revisited 

Taking into account that 8 f2 (E) = Y^=o ^« {E), and using (2.6), for the single string 
density of states when E > A, Bra ndenberger and Vafa obtained 

6 Only for the Maxwell-Boltzmann statistics fi„ can be understood as the density of states for a 
gas with n strings [2]. 

7 This means that we have used an exact form for the coefficients giving the degeneration number 
at each mass level of the superstring and the sum over winding and momenta has been obtained by 
integration. 

8 For n — we have that f2o (E) corresponds to the vacuum state, whose density of states is 
Dirac's delta. 




(4.1) 
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ay z (56, t) 




Figure 1: 002 (56, t) computed numerically. The continuous line represents 0J2 as given by 
using Q'l only (a' = 1). 



where: 



1 /" , fo ^ cq SI -i / /" , sinx \ 
— / clae^o d *— ^- cos / dx— a = 0.56 ± 0.01 



-00 \o 



00 



1 / , f« j^ cosx-i I f , sinx \ 
6 = — / «dae J » Ax ^r~ sin / dx a J = -0.29 ± 0.01 

-00 \o / 

The values of a and b were calculated numerically. Note that, although the exponent 
gives, when a>l,a factor — logo, the second integral is not convergent because of 
the oscillatory behavior of the integrand, and so it has to be regulated and they did it 
by using an exponentially decaying factor. The sign of b produces a positive specific 
heat, and the Hagedorn temperature is then a maximum one in the micro canonical 
ensemble. We think, nevertheless, that this scenario needs a revision: there is nothing 
a priori unphysical in negative specific heats in the microcanonical ensemble, and 
anyhow, we are not sure that there could be a positive specific heat phase in the 
microcanonical ensemble for the gas of strings with windings and momenta. 

In fact, it is easy to see how a contradiction appears, when b = —0.29, in the fol- 
lowing way: using (4.2) and the fact that Z = £ {Q (E)} the following expression 
for the partition function can be obtained 9 : 



9 The multiple string density of states must have a cutoff that indicates the range of validity of 
the asymptotic approximation and that would be, in general, different from the single string cutoff. 
We are using <fi as the cutoff for the multiple string density of states. 
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Z h (P) = A(f3 _ ^ + ^ [0,0 {(3 - (3 H )] (4.3) 

whereas from the single density of states it is possible to arrive also to an expression 
for the multiple string partition function via the equality — f3F ((3) = C {VL\ (E)} 

Z\(3)= e-W) = — + + O [(/? - (3 H f] (4.4) 

Comparing (4.3) and (4.4), one immediately gets that 10 a = exp (— 7) = 0.5614, 
which perfectly agrees with the numerical value given in [1]. But it is a very notorious 
fact that in (4.4) no term is found analogous to the logarithmically divergent one 
hidden in the incomplete gamma function of (4.3). 

Furthermore, we have found an alternative calculation where the coefficient b is 
actually zero and then, only one divergent term is present in both equations (4.3) 
and (4.4). As it is written in (2.6), the high energy dominant term in the single 
string density of states has a dimensionless factor which equals unity in the type II 
and the heterotic string (this factor could be volume dependent when considering 
open strings and branes). However, it is very useful to introduce in Qi (E) a factor c 
that could be thought just as a regulator (whereas we are going to give it a physical 
meaning at the end of this subsection). This change adds a factor c n to f2 n in (4.1). 
Only at the end the limit c approaching to one will be taken. Furthermore, it is 
possible to work out the values of a and b analytically if a simple change is made in 
(4.1): the key ingredient is noting that the upper limit in the integrals can be taken 
to infinity since the Dirac delta function ensures that no value greater than E will 
contribute to them. This will render the final integrals much easier to perform; then, 
we will have 



E +00 +00 

dEi — / dEi n h (E) = ^ / dae- ia e ^ * (4.5) 

A A -00 

The calculation is now analogous to the one made in [1], taking the first terms 
in the series expansion in A/E. 



P„E [ roo-sg(a) dx j r oo-sg(a) 

I i — irv J a A x J oA x 

e ~e~ e "b" 



—00 
+00 



2ttE 

—00 



J da e- Ci (^) cos (a + c • si (^j ) (4.6) 



10 This expression for a was roughly deduced in [11], although no connection with the numerical 
value of [1] was made there. 



- 10 - 



Where sg (a) stands for the sign function. It is important to note that the Ci and the 
si functions only coincide with the standard cosine integral and sine integral functions 
for positive values of a. Ci (x) is a real, even function of x whereas si (x) is a real, odd 
function of its argument with a discontinuity at x = and si (0) = 0. This is easy to 
understand as a consequence of the integral upper limits including the sg (a) term in 
the exponential factors in (4.6). The integrand is then even and one could perform a 
series expansion in terms of aA/E. Approximating Ci (^) = 7 + log (^) + O (^) 
and si (^) = =f + ^ + O {^-f (a > 0), we arrive at 



00 



ohms _ 7 e feE /£V (da / tt aA\ 
Q n (E) = e 7 — — — / — cos [a - c- + c— (4.7) 
K J nE \A / a c V 2 E 1 







Using that cos (a — b) = cos a cos b + sin a sin b and taking the lowest order in 
aA/E one finally gets 



0^(£) = — — e 
v 7 A c 



<c) + 6(c) A + O (± 



(4.8) 



00 

~ n r da ( cix 



a (c) = / — cos ( cm = - . . if < c < 1 

w n a c V 2 / T (c) 





00 



, / n ce n f da . (cix \ ce 01 

b (c) = / sin a) = — if 1 < c < 2 

K ' 71 J a c ~ l V 2 / r(c- 1) 


The integrals converge only for the indicated values of 11 c. Clearly, in the c — > 1 
limit, the result for a is fully compatible with the numerical value given in [1]; the 
problem is that it is very easy to see how in this limit b goes to zero. As a matter of 
fact this method can be generalized and more terms of the form d n (A/E) n (n G N) 
can be computed, giving all of them zero for d n in the c — > 1 limit. 

Now, it is straightforward to see how the contradiction between the equations 
analogous to (4.3) and (4.4), that now depend on c, has disappeared. Using fix (E) = 
ce^ HE /E we have that 



Z h {(3)= e~ c/3Fh ^ » *~ C \ - c + C6 " C7 - , + O {{(3 - l3 H f~ c ) . (4.9) 

And making directly the Laplace transform of (4.8), the same expression for Z (P) = 
C {fl (E)} can be found for c > 1. When c = 1 an annoying constant term appears 
preventing us to fix more than the only divergent term. 

11 One could be tempted to put directly, in the expression for a (c), that, when c = 1, cos (| — a) = 
sin a; but this would produce a wrong result since we would be forgetting that si (0) = 0. 
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Once the value c = 1 is taken, the expression for the density of states of the 
string gas in a finite size container is given by 

n h (E)=e~^^- £>A. (4.10) 

From this we can calculate both the entropy of the system and its temperature. 
The fundamental thermodynamic relationship giving the entropy as a function of the 
energy now looks like 

S h W=j3 H E + ]ag(jp). (4.11) 

Comparing it with (2.10), the analogous expression in the fixed temperature, case 
we see how both ensembles seem to be inequivalent as pointed in (2.8). With this 
entropy we will also have that temperature is fixed to Hagedorn's, and we would 
have to conclude that Cy (E) would be infinite. 

The constant c has been introduced as a mere way of doing analytical continua- 
tion of ill-defined expressions, but we can give it a physical interpretation. Lets look 
at the expression defining Q (E) once c is introduced 

oo 

n(E) = J2c n ^n(E) (4.12) 

?1=0 

We see that c is really acting as the fugacity of the system, so that working with 
a generic value of c and then performing the c — > 1 limit is exactly the same as 
working with a generic non null chemical potential and then taking it to zero. This 
interpretation also lets us know that we have not been working in the microcanonical 
ensemble, but in the "enthalpic" one [2] for which energy and the chemical potential 
are given. This way Q (E) now depends on c and becomes Q (E, c). 



4.2 Fluctuations for the fixed energy description 

It is now clear that the density of states of the string gas is given by 

n h (E) = ^-e(E-(j>)^ E (4.i3) 

The multi-string energy cutoff <p cannot be completely determined without matching 
the high energy regime with the low energy phase because, imposing that f3F h (j3) 
must be obtained, <p would appear as a factor of (/3 — (3h)° that is a regular term 
being F l (f3) also regular at f3 H . The matching can be done, but we are not going to 
dwell further on this point. 
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An immediate consequence of (4.11) is that the high energy microcanonical spe- 
cific heat, Cy (E), is divergent contrary to what, by using [1], has been assumed as 
true for more than fifteen years 12 . 

Once Q (E, c) for the enthalpic ensemble is obtained, it is straightforward to 
compute the number of strings and its fluctuations as 



N(E) = cd c \o g n(E,c)\ c=1 ^log 



AN 2 (E) ^ {cd c ) 2 log fl{E,c) 



c=l 



N(E) 



N(E) 



N(E) 



with 



Q h (E,c) = d(E-<P) 



-cy 



A c r c 



l + (c-l)0(-)+... 



A 



that has been obtained by Laplace inversion of e cq ^\ 



(4.14) 



4.3 Other refinements 

In the preceding sections any dependence on R has been lost as a result of being in 
a physical situation in which sums can be well approximated by integrals. Now, we 
would like to add the effects of introducing Euler-Maclaurin corrections [4] in the 
integrals that represent sums over windings and momenta. This can also be done by 
means of a Poisson resummation. As a result, there appear more singular points in 
—(5F{j3) whose location depends on R [11]. 



-pF h (p)=J2 r 



with 



m,n,j 



o, 



2tt (p*-fa Hij (R) 2 ) 



^m,n, 3 (j + 1/2) 2 



R'' 



a' 



2 -m — — n' 



a' 



R- 



(4.15) 



(4.16) 



Assuming R > \fa' the singularity nearest to (3 H is p 1 = f3 H — rj, rj ps a'(3 H /(4R 2 ). 
Considering also the term depending on (5\ we get a corrected expression for the 
density of states 



12 In the next subsection and the conclusions, we will treat and discuss the relevance that the 
radius corrections can have in relation to this thermodynamical statement. 
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«> = (l - (g) " ,4,7) 

Where 77-E 1 ^> 1 has been assumed. This expression depends on R through 7] and the 
R corrections are actually non extensive. If taken into account, the effect of these 
corrections on j3 would be 

[3 (E) = fa + (vE) 17 (4.18) 

The second term on the right hand side of this expression has already been stud- 
ied in the literature [11], [12] and has been frequently claimed as being the cause 
of having a positive specific heat. It is easy to see that, after all, having enough 
energy for a given radius, the gas can always be in the regime in which there is no 
dependence on the volume and this is the regime of the Brandenberger-Vafa scenario, 
for which the specific heat diverges. The radius corrections are exponentially sup- 
pressed and, for our system, are not different from the finite volume corrections that 
are dropped when the thermodynamic limit is taken over the ideal gas of particles 
(see, for example, [2]). 

5. Conclusions 

As a first important result, it is crucial to remark that no published calculation 
has found the —0.29 e l3nE /E term but the one presented by Brandenberger and 
Vafa in [1]. This is the term that is needed to state, as these authors do, that 
the microcanonical specific heat is positive in the physical situation in which energy 
density is so high that the density of states does not depend on the volume. As far 
as we know, what any other calculation actually gets is that, for the same physical 
situation, the specific heat is divergent because a null b coefficient is found (see 
eq. (4.2)). In reference [11], it is explicitly admitted that b = —0.29 is not found, 
but the authors do not face the important question that the contradiction between 
their calculation and that of Brandenberger and Vafa rises. It seems that this is 
so because they consider they are using what they call a different "formalism". In 
our work, we clearly show that, using the same technique Brandenberger and Vafa 
used and a physically meaningful regularization 13 , the b coefficient vanishes (and also 
any other high energy correction). This is the first result of our work and, in our 
opinion, critically depends on understanding that the " micro" ensemble is really the 
"enthalpic" one, namely, a fixed energy and fixed chemical potential ensemble. 

13 Brandenberger and Vafa stay that they use an exponential regulator for a numerical computa- 
tion. 
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We have used the ^-representation of the Helmholtz free energy to finally con- 
clude that the behavior of the free energy around (3h coincides with what is gotten 
from the F representation. This might be expected but it is not a trivial fact 14 , 
because both representations do not coincide. That the S and F representations 
are not equal seems to be clearly commented for the first time in [13]. The re- 
lationship between both representations is carefully studied for the family of the 
heterotic strings in [14] (see also [15]). In those works it is clearly established that 
the F-representation does not provide an analytical function for complex (3. In other 
words, it is false that the F-representation can be seen as providing the analytical 
continuation of the ^-representation for heterotic strings. Then, it is not rigorous 
to say that the F-representation of the free energy can be used to get the density 
of states by inverse Laplace transformation of the corresponding partition function. 
What one can only do is to continue the ^-representation. When one uses the F- 
representation to get the behavior around (3h one is really using it in the interval 
(Ph,+oo) where it coincides with the ^-representation. Now, we are also providing 
a very concrete and explicit example of to what extent the S and F representations 
of the free energy are equal. 

From other point of view, taking into account that the F-representation gives 
the free energy as computed for the compactification of the Euclidean time on a 
circle of length f3 (including string windings along it), in [16] there appears a proof 
for the noncritical c = 1 string of the recently emphasized fact that, for strings, the 
free energy cannot be computed, at any temperature, by the compactification of time 
[17]. 

Another important point is to what extent the R corrections presented in sub- 
section (4.3) can be taken as showing that the system really has a positive specific 
heat. This seems to be the belief as expressed in [11] and [12]. In our opinion, the 
problem is that these nonextensive corrections are exponentially suppressed with the 
value of the energy and are then of no thermodynamical relevance for the system 
in the thermodynamical limit in which the energy density is very high and the sum 
over windings and momenta can be replaced by an integral. Those corrections are as 
the finite volume l/W corrections for the gas of free particles; they are irrelevant 
in the thermodynamical limit in which the sum over momenta can be replaced by an 
integral. 

In the case we treat the gas in the fixed temperature (canonical) description, the 
big fluctuations might justify the exclusion of the nonextensive term that renders the 
canonical specific heat positive. In any case, those fluctuations would just make the 

14 In the previous versions of our work we thought on the contrary because we found only the 
leading order contribution in (2.2). In fact, a seemed very concerned referee, after reading the first 
version of our work, was fully convinced of the difference between the F and S representations for 
this calculation and used it to reject the work without any hesitation because he/she liked more 
the F-representation. 
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canonical equilibrium description physically unusable. 

From a cosmological point of view, one could think that there would not be any 
relevant cosmological implication from our results because, after all, the equation of 
state would still be the same, corresponding to pressureless dust matter (P = 0). 
However, things are more intricate because a divergent micro canonical specific heat 
can be an indication of a phase transition. In our case, what we have done in 
the microcanonical (really enthalpic) treatment is a description of how the volume 
of phase space in the N-body problem changes when energy, which is a conserved 
quantity, increases [18]. We have found that the system behaves very differently from 
the grand canonical ensemble in which Tjj would be a maximum temperature and 
the specific heat, as a function of temperature, would be positive. 

What is clear is that a divergent microcanonical specific heat cannot be used ab 
initio as a criterion to drop as unphysical our gas of closed strings at finite size. 
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A. The UV limit in the S-representation partition function 

This appendix is devoted to explain how (2.2) is obtained and to explicitly show that 
(2.4) holds. 

First of all, the integral computing the sum over windings and momenta can be 
performed to give 

J dfd£e-^ f2+p )e 27ri? * T1 = \t\- 9 (A.l) 

R s xR 9 

Next, it has to be noticed that, when r — > 0, it holds that 

|^ 2 (0,-l/(2r))r 16 ^2- 16 e 2 ^ 2 /l T l 2 

because all the other terms are finite when r — > 0. One has then to perform the 
integral over T\ as providing a function of r 2 given by e 27r / T2 times a series expansion 
in powers of r 2 as it appears on the right hand side of (2.2). Namely, the left hand 
side of (2.2) (we called it I(t 2 )) is now given by 



2 8 r 2 



+1/2 

3 2 " /T2 J dn 

-1/2 



-1/2 



T ■> 

"2 



(A.2) 
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It has been written in a way prepared to be rewritten in terms of a function /(r 2 ) as 
J(r 2 ) = 2- 8 r 2 1/2 e 2 ^ T2 J(r 2 ) where 

J(r 2 ) = J dx e~ 2 ™ 2 ettiT^-l 1 -^ 2 )] (1 + r 2 a; 2 )" 1/2 (A.3) 

For which the change of variables T\ — x t^ 2 has been used. Now, for the product of 
the last two factors in the integrand, the following series expansion can be written 



fi b^-t 1 -^ 2 )] (1 + T2 x 2 y 1/2 = e 2 - 4 W(i+^ 2 ) (1 +T2 x 2 y 1/2 

a (27r) 6 r(6 + a + 
6! a! r(6 + l/2) 



^ ^ (-l)'( 2 7r)»r(6 + a + l/2) b+a fc+2a 



6=0 a=0 

4fe+2a 



Next we are able to perform the integral over x of the term x + obtaining 



KT 1 



/ 



^-2^46+20 = F + a + V2) ~ F + a + 1/2 ' ^ ~> /(4T23)) 



( 27r )2fe+a+l/2 

-KT 1 

^ (2n)- 2h - a - 1 / 2 T(2b + a + l/2) (A.5) 

where the last approximation results from the exponential suppression of the con- 
tribution coming from the incomplete Gamma function when its argument gets big 
(what happens here when r 2 — > 0). 

Now one of the sums in the resulting double sum to get I can be calculated (!!) 
to finally give 



J(r 2 ) = (27T)- 6 - 1 / 2 r (b + 1/2) r 2 b (A.6) 
6=0 

We are then able to get /(t 2 ) and, in particular, the coefficients a« as 

r (i + 1/2) 

ai = (A ' 7) 

The next step is to use /(r 2 ) as written in (2.2) with the already known aj coefficients 
to compute F h ((3) finding then the b n factors in (2.3). This is easily done to give 
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The final computation is that of the series generated by the b n coefficients as it 
appears on the left hand side of (2.4), namely 



+00 



n=0 



(A.9) 



Taking into account that (3 2 - (3% = 2(3 H ((3 - (3 H ) (l + ), R{(3) can be written 
as a double sum 



+00 



m = E 



n=0 



1 n 

77^ E 



-iy2 2q ~ n T (g + 1/2) 



;r ^5 (n-g)!(2g-n) 



(0 - /3^) n 



)9h 



n+l 



(A.10) 



The term between square brackets can be computed to give exactly (— l) n . So, we 
have finally showed that 



R{(3) = 1/(3 (A.11) 
by showing that it is given by the power series expansion of 1/(3 around (3h- 
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